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Band insulators appear in a crystalline system only when the filling - the number of electrons per unit cell 
and spin projection - is an integer. At fractional filling, an insulating phase that preserves all symmetries is a 
Mott insulator, i.e. it is either gapless or, if gapped, displays fractionalized excitations and topological order. 
We raise the inverse question - at an integer filling is a band insulator always possible? Here we show that 
lattice symmetries may forbid a band insulator even at certain integer fillings, if the crystal is non-symmorphic 
- a property shared by a majority of three-dimensional crystal structures. In these cases, one may infer the 
existence of topological order if the ground state is gapped and fully symmetric. This is demonstrated using 
a non-perturbative flux threading argument, which has immediate applications to quantum spin systems and 
bosonic insulators in addition to electronic band structures in the absence of spin-orbit interactions. 



Fundamental to the study of crystalline materials is the ex- 
istence of band insulators, and the fact that they only occur 
when the electronic filling is commensurate with the lattice. 
Defining the filling as one half the number of electrons in 
a unit cell - which accounts for spin degeneracy in the ab- 
sence of spin-orbit coupling - band insulators occur only if 
the filling is an integer. At fractional filling a material must be 
metallic within the free electron approach. However, interac- 
tions could lead to a distinct insulating phase. For fermionic 
systems, this would have to be a Mott insulator, strictly de- 
fined as an insulator that breaks no symmetries but is distinct, 
i.e. separated by a phase transition, from a band insulator. 
In fact, a far stronger result is true: even interactions can- 
not cause a system at fractional filling, whether of bosons or 
fermions, to enter a trivial insulating state. There are only two 
choices if all symmetries are to be preserved: the system must 
remain gapless, or develop topological order, which we take 
here to mean a ground state degeneracy that is not associated 
to a broken symmetry, as exemplified by fractional quantum 
Hall states and gapped quantum spin liquids. These gapped 
states feature long range quantum entanglement, leading to 
excitations with anyonic statistics and ground state degener- 
acy that is sensitive to the spatial topology. A concrete ex- 
ample is the the one-dimensional (d = 1) Hubbard model at 
half filling, with repulsive interactions whose ground state is 
insulating. The spin excitations in the insulator form a gap- 
less Luttinger liquid, as expected of a fully symmetric state in 
d = 1, where a gapped state with topological order is forbid- 
den. This result was demonstrated for arbitrary spatial dimen- 
sion by Hastings 1 2 and Oshikawaj^by extending the d = 1 
Lieb-Schultz-MattiP'theorem to d > 1 using a beautiful 'flux- 
threading' argument.^ 

However, surprisingly little is known about the inverse 
question: given an integer filling, is it always possible to find 
a 'trivial' insulating state, like a band insulator, that is not re- 
quired to have these exotic properties? Here we show that for 
many crystals, trivial insulators can be forbidden even though 
the filling is an integer. For bosons or fermions, if the insu- 
lating states preserve symmetry, they must either be topolog- 



ical^ ordered or host gapless excitations. We demonstrate a 
simple crystallographic criterion for this to hold: namely, that 
the space group of the crystal is non-symmorphic. A crys- 
tal is non-symmorphic if there is no choice of origin about 
which all its symmetries can be decomposed into a product 
of a lattice translation and a point group symmetry element. 
Our key result is succinctly stated: There are integer fillings 
at which insulators with unique ground states are impossible 
if the space group is non-symmorphic. For such crystals we 
show that there is a minimal integer filling, strictly greater 
than unity, at which trivial insulators are allowed. We call this 
the non-symmorphic rank, denoted by S > 1. Trivial insula- 
tors exist only at fillings which are integer multiples of S. 

We may rationalize this unexpected result by observing that 
in non-symmorphic crystals, a fractional lattice translation 
acts in concert with another transformation to leave the crys- 
tal invariant. This, loosely speaking, renders an integer filling 
equivalent to a fractional one, thus forbidding a trivial insula- 
tor. Similar conclusions can be drawn for any system where it 
is possible to identify a conserved charge. These include spin 
systems with a conserved spin component, and lattice bosons 
whose number is conserved. Note, 157 of the 230 distinct 
space groups in three dimensions are non-symmorphic and 
have S > 1. These include the ubiquitous hexagonal close- 
packed (hep) structure and the diamond lattice space group 
relevant to many semiconductors and pyrochlore materials. 
Table [I] has further examples. The implications therefore are 
manifold, a few applications will be discussed below. 

Before we move to the main results, it is worth making two 
observations. First, we make no assumptions beyond lattice 
symmetry. In particular we do not limit ourselves to tight- 
binding models where additional restrictions may emerged 
Indeed, the construction of trivial interacting insulating states 
in those cases relies precisely on the fact that they do not vi- 
olate the conditions we establish here. Second, when we dis- 
cuss electronic band structures, we assume that the electron's 
spin is purely a spectator, i.e. there is no spin-orbit coupling. 
This allows us to address each spin species independently, as 
though the electrons were 'spinless'. Moreover, we will con- 
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nect our key results to extinguished Bragg peaks, which as- 
sumes a 'scalar' coupling between electrons and the crystal 
lattice. Incorporating electronic spin-orbit coupling is an im- 
portant extension which we leave to future work. 



I. APPLICATIONS 



(i) Band Theory.- The fact that band insulators are forbidden 
unless the filling is a multiple of the non-symmorphic rank 
S, strongly constrains the structure of the bulk energy disper- 
sion: it is impossible to 'detach' a set of fewer than S bands 
so that they touch no other bands, without breaking the crys- 
tal symmetry. Non-symmorphic space groups often also de- 
scribe the symmetry of photonic crystals, 9 10 hence their pho- 
tonic band structures obey similar constraints. For instance, 
the hep structure has S = 2; we show a tight-binding band 
structure for the hep lattice in Fig.|2]where the enforced con- 
tacts are explicit. Although the subject of band touchings in 
crystals has a long historypE^l an( j the ubiquity of such de- 
generacies in non-symmorphic crystals has been notedP^EI 
the connection to a minimum filling for band insulating be- 
havior has not been made previously. More importantly, these 
prior results apply only to noninteracting systems, in contrast 
to the non-perturbative approach taken here, which allows us 
to determine the nature of the interacting insulating ground 
state. 

(ii) Spin Systems and Bosonic Insulators.- A parallel set of 
conclusions can be drawn for spin systems, where the filling 
is related to the total magnetization, and an insulator now cor- 
responds to a phase with a spin gap. We demand at least a 
U(l) spin rotation symmetry, although the conclusions also 
apply to a larger symmetry such as SU(2), so long as it con- 
tains a U(l) subgroup. The analog of the band insulator is a 
trivial paramagnet, that has gapped excitations, and has nei- 
ther conventional nor topological order. We conclude that 
such a trivial paramagnet is disallowed in an SU(2) symmet- 
ric spin- 1/2 system on the diamond lattice, which is a com- 
mon non-symmorphic lattice with 6> = 2. In contrast, the 
pyrochlore lattice has the same space group but twice as many 
sites per unit cell as diamond, and a trivial quantum param- 
agnet is not forbidden by our arguments. Another application 
is to magnetization plateaus in an applied Zeeman field. For 
example, half magnetization plateaus of spin- 1/2 moments 
on the two-dimensional non-symmorphic Shastry-Sutherland 
lattice 15 (SSL) cannot be trivial paramagnets. Applications to 
SrCu 2 (B0 3 ) 2 (SCBO), a material realizing the SSL where a 
half magnetization plateau is observed] 16 ! 17 ! will be discussed 
below. Yet another application is to bosonic Mott insulators. 
Our arguments demonstrate that on non-symmorphic lattices, 
Mott insulators at fillings that are not a multiple of S must 
be topologically ordered, if they are gapped and respect all 
symmetries. 



II. FLUX THREADING ARGUMENT AND 
NON-SYMMORPHIC RANK 

In the balance of this paper we outline our argument, and 
demonstrate its use in specific examples by substantiating the 
claims made above. For clarity and to fix notation we digress 
briefly to review some relevant crystallography. We will con- 
sider crystalline systems with a given space group, Q. This 
has two ingredients: (i) the subgroup of translations T, gen- 
erated by the set of primitive translations i a . : r — > r + a i7 
i = 1,2,3, where the generate a specific Bravais lattice, 
and (ii) a point group, V consisting of discrete rotations, in- 
versions, and reflections. By combining the 14 Bravais lattices 
with the 32 crystallographic point groups, we obtain the 230 
space groups. Of these, 73 are symmorphic: there exists at 
least one point that is left invariant by all the symmetries, up 
to translations by a lattice vector. The remaining 157 space 
groups for which no such point exists, are non-symmorphic. 
Intuitively, a non-symmorphic space group contains one or 
more essential non-symmorphic operations ('glide mirrors' or 
'screw rotations') that combine a point-group operation with 
a fractional lattice translation, that cannot by any change of 
origin be rewritten as the product of a point-group operation 
and an ordinary translation. The latter caveat is important to 
distinguish these from 'trivial' glides/screws which can occur 
even in symmorphic crystals. 1 - There are however two 'ex- 
ceptional' space groups, which are non-symmorphic despite 
the absence of essential glide or screw operations™ 

In the remainder, we focus on the 155 non-symmorphic 
groups which have glide mirrors or screw rotations. Under 
a non-symmorphic symmetry G consisting of a point group 
operation g and non-lattice translation t, the Fourier com- 
ponents of a scalar field (such as the density) transform as 
rife — > ngke ll ~' k . If gk = k and r • k is not an integer multiple 
of 2ir, rik — 0: the associated Fourier component vanishes. 
An essential glide or screw always has an infinite set of such 
reciprocal lattice vectors k. Thus non-symmorphicity (barring 
two exceptions) has a dramatic experimental manifestation: 
there are systematic absences in the diffraction pattern where 
the (scalar) Bragg intensity vanishes. 

We assume that the systems we study are described by a 
Hamiltonian H that preserves all the symmetries of Q and in 
addition that there is a conserved U(l) charge, Q (assumed 
to take a fixed integer eigenvalue and thus replaceable by a 
c-number throughout) with [H, Q] — 0. We make no as- 
sumptions as to the origins of the conserved charge, so for 
instance the systems we consider could be built out of (i) spin- 
less fermions or bosons, where Q is just the conserved parti- 
cle number; (ii) spinful fermions with SU(2) spin symmetry, 
in which case Q is one-half the total fermion number (since 
the two spin components may be treated independently); or 
(iii) lattice spins with (at least) U(l) spin rotational invari- 
ance, in which case we may take the charge on lattice site r 
to be S + where S, S* are its total spin and magnetization, 
and define Q accordingly. Considering a finite system with 
Af c = N 3 unit cells, we may then define the filling to be the 
charge per unit cell, i.e. v — Q/Af c , which will be held fixed 
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(a) (b) 

FIG. 1. Flux Threading on a Non-symmorphic Lattice, (a) A 2D 
non-symmorphic lattice, with p4g space group. An essential glide 
symmetry is shown, consisting of a mirror (02) and half-translation 
(T ai /2)- (b) Flux threading changes ground state momentum from 
Pi to Pf. At odd integer filling, Pf is an extinguished reciprocal 
lattice vector, denoted by empty circles, implying a distinct ground 
state. At even filling, Pf lands on an allowed reciprocal lattice vec- 
tor, allowing for a unique ground state. 



in the thermodynamic limit. We impose periodic boundary 
conditions that identify r and r + iVa*. 

Our argument is based upon the response of the system to 
the insertion of a gauge flux that couples minimally to the con- 
served charge Q. We first give a heuristic argument, and then 
a more formal one. The strategy is as follows: we wish to 
show that if the system is an insulator, the ground state cannot 
be unique in the thermodynamic limit. Hence the system can- 
not be a trivial insulator that respects all symmetries, which 
should have a unique ground state. A degenerate ground state 
implies either (i) a broken symmetry or (ii) a gapless phase or 
(iii) topological order. That is, it is topologically ordered if it 
is gapped and breaks no symmetries. To show this, we begin 
with a ground state \ty) and thread a flux quantum through a 
periodic direction, which, by gauge invariance, returns us to 
the original Hamiltonian. This procedure produces an eigen- 
state |$). Earlier worlsPSliE] has argued that for an insu- 
lator, \^f) must be a iow energy' state, i.e. its energy ap- 
proaches that of the ground state in the thermodynamic limit. 
Although rigorous energy bounds can only be given for a dif- 
ferent - but gauge-equivalent - flux insertion, for pedagogical 
reasons we keep a simpler choice with the understanding that 
the arguments of Ref. 2 can be applied, mutatis mutandis, to 
non-symmorphic symmetries. The key step is to show that 
|\&) is distinct from |^>), which would then establish ground 
state degeneracy. In the case of fractional filling, these states 
differ in crystal momentum ! 1 l 2 l 2 ° l 2 1 1 For integer filling, crys- 
tal momentum fails to differentiate between them. However, 
we show that on non-symmorphic lattices, one can still dis- 
tinguish these states using the quantum numbers of the non- 
symmorphic operations. 

Consider threading flux through the system by introducing 
a vector potential A = k /N , where k is a reciprocal lattice 



vector; this is the most general vector potential that is 'pure 
gauge' so that e l $c A dr — 1 for any loop C that threads the 
system. The change in energy upon inserting flux is propor- 
tional to the total current and is thus bounded in an insula- 
tor, where this procedure produces a state degenerate with the 
ground state in the thermodynamic limit. If the flux threading 
changes the quantum numbers, the final state is distinct from 
the initial one and thus the ground state is degenerate. We use 
units where H = e = 1, in which the flux quantum 0o = 2ir. 
In the case of spinful fermions, we consider coupling to a sin- 
gle spin species, justified by the fact that spin is assumed to 
be a spectator. 

In order to determine the change in quantum numbers of 
the symmetry operators, it is useful to first ignore the lat- 
tice potential, in which case it is straightforward to identify 
the change of the center-of-mass momentum by computing 
the force imparted to each charge as the flux is adiabatically 
switched on. Faraday's law gives F{ = J A ■ dr, so that 

a r n 
AP = ^ / Fidt = ±k = N 2 vk (1) 

2 — 1 

We now reintroduce the lattice. If AP is not in the reciprocal 
lattice, then it is observable even within the reduced symmetry 
of the crystal, and we have succeeded in producing a distinct 
degenerate ground state. For a fractional filling v = p/q, 
it is clear that this is the case so long as we choose N rel- 
atively prime to q. In other words, the state following flux 
insertion has a distinct crystal momentum, which means that 
the quantum number associated with translational symmetry 
has changed. This is the essence of the Hastings-Oshikawa 
argument. At integer filling, it is clear that no choice of N 
allows us to distinguish the initial and final states on the basis 
of lattice translational symmetry. 

However, in a non-symmorphic crystal, if we can choose 
AP to correspond to an extinguished Bragg peak, the ini- 
tial and final states will have different quantum numbers for 
the non-symmorphic operation responsible for that extinction. 
Such a choice is always possible at unit filling, implying that 
an insulator with a unique ground state is forbidden. Intu- 
itively, the absence of Bragg peaks at special crystal momenta 
indicates that they can serve as a good quantum number to dis- 
tinguish initial and final states, in the presence of a crystalline 
potential. In general, as long as v is not a multiple of the 
non-symmorphic rank S, a unique ground state is forbiddenfSl 
This result may also be understood by the exercise of trying 
to isolate a single Bloch band in a non-symmorphic crystal 
without breaking symmetry. Were such a band to exist, its 
exponentially localized Wannier orbitals would define centers 
of electronic charge in each unit cell. One now encounters an 
obstruction: it is impossible to define a charge center that is 
invariant, modulo translations, under all crystal symmetries. 

We now bolster this intuitive picture with a more formal 
argument. Consider a crystal which has a non-symmorphic 
space group Q, which contains a non-symmorphic operation 
G. This comprises a point group operation g followed by a 
fractional translation r in a direction left invariant by g, i.e. 

G :r -t gr + r (2) 
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where r is not a lattice vector, and gT = r. We begin with 
a ground state \ty), and assume it is an eigenstate of all the 
crystal symmetries, including G, i.e. 

G\%)=e ie \$>) (3) 

We consider the smallest reciprocal lattice vector k left in- 
variant by g, so that gk = k and k generates the invari- 
ant sublattice along k. We now thread flux by introduc- 
ing a vector potential A — k/N as before, in the process 
of which \^>) evolves to a state \^f') that is degenerate with 
it. To compare to \*ff), we must return to the original 
gauge, which can be accomplished by the unitary transforma- 
tion I*') £7 fc |*') = |#), where 

U k =expS.^ J d d rk-rp(r)\ (4) 

removes the inserted flux, and p(r) is the density operator cor- 
responding to the conserved charge Q. Since A is left invari- 
ant by G, threading flux does not alter G eigenvalues, so 
and \ty') have the same quantum number under G; however, 
on acting with U k , the eigenvalue changes, as can be com- 
puted from the equation: 

G- 1 U k G = U k e 2 ^^ Q / N (5) 

where we have defined the phase factor <& s (fe) = r-fc/27r,and 
Q = vN 3 is the total charge. It may be readily verified that 
since gk = k, & g (k) is unchanged by a shift in real-space ori- 
gin. For a non-symmorphic symmetry operation G, this phase 
<& g (k) must be a fraction. This follows since r is a fractional 
translation. (If a lattice translation had the same projection 
onto k as r, this would yield an integer phase factorial How- 
ever, this would render the screw/glide removable i.e. reduced 
to point group element x translation by change of origin.) 
Thus, for G non-symmorphic, & g (k) = p/Sq, with p, Sg 
relatively prime. From ^ we conclude that |^) and |^) have 
distinct G eigenvalues whenever $ g (k)Q/N — pN 2 v/Sc is 
a fraction. Since we may always choose N relatively prime 
to the Sg, the result of flux insertion is a state distinguished 
from the original state by its G eigenvalue, unless the filling 
is a multiple of Sg ■ For a glide Sg = 2, while for a screw Sg 
is the number of times it must be applied before it becomes 
removable.^! 

Non-Symmorphic Rank - Recall that trivial insulators are 
only allowed at filling that are multiples of an integer S, which 
we call the non-symmorphic rank. Each non-symmorphic op- 
eration G is individually associated with an integer Sg which 
leads to a degeneracy unless it divides the filling v. Hence the 
non-symmorphic rank S is divisible by the least common mul- 
tiple of the Sq- A tighter bound on the non-symmorphic rank, 
although subtle to prove^is easy to state: it is the smallest in- 
teger S such that some point is invariant under all the elements 
of Q up to 1/5 times Bravais lattice vectors. 

Examples. -The statements made earlier about specific ex- 
amples follow immediately from an examination of the 
space group symmetry!^! (i) The hep crystal (space group 




(a) (b) 



FIG. 2. Band Structure on hep lattices, (a) hep crystal showing 
screw rotation comprised of sixfold rotation (R^ ) followed by half- 

c-axis translation (T c / 2 ) out of the page; green arrows show transfor- 
mation of some sites, (b) (bottom) Tight-binding band structure (blue 
solid lines) showing symmetry-enforced contact^ between the pair 
of bands, leading to a minimum filling of S = 2 to achieve a band in- 
sulator, (top) On breaking the screw symmetry, gaps open and allow 
the bands (red dot-dashed lines) to separate. 




FIG. 3. Other Examples of Non-symmorphic Lattices, (a) The dia- 
mond lattice space group FdSm has a screw rotation comprised of 
fourfold rotation (R n /2) followed by a quarter-translation along z 
(Tg/i). (b) The non-symmorphic Shastry-Sutherland lattice has the 
same p4g space group as the lattice of Fig.[TJ 



PQz/mmc) possesses a sixfold screw rotation about its c- 
axis, two applications of which result in a pure translation; 
as the only other non-symmorphic operation is a glide also of 
rank 2, we conclude S = 2 . A model band structure, with 
enforced contacts between pairs of bands is shown (Fig. [2}; 
breaking the symmetry explicitly lifts these degeneracies, de- 
picted by the broken red lines, (ii) Diamond and pyrochlore 
(space group Fd3m) possess a four-fold screw axis, but ap- 
plying this twice yields a removable twofold screw (Fig.[3](a)). 
Once again, all remaining non-symmorphic operations have 
rank 2, so that 5 = 2 for both these crystals, (iii) In a field, 
SCBO exhibits a magnetization plateau at half the saturation 
va i ue ji6liZl J \s this corresponds to v = 1, unlike other plateaus 
which are at fractional filling j 22 " 25 ' Hastings-Oshikawa argu- 
ments cannot be applied to it. However, the two-dimensional 



5 



d 


Name 


Examples 


Space Group 


5 


2 


Shas try-Sutherland 


SrCu 2 (B0 3 )2 


pig 


2 


3 


hep 


Be, Mg, Zn 


PGa/mmc 


2 


3 


diamond 


C, Si 


Fd3m 


2 


3 


pyrochlore 


Dy 2 Ti 2 07 (spin ice) 


Fd3m 


2 


3 


hyperkagome 


Na 4 Ir 3 08 


P4i32 


2 


3 




a-Si0 2 , Ge0 2 


P3i21 


3 


3 




CrSi 2 


P6 2 22 


3 


3 




Pr 2 Si 2 7 ,La 2 Si 2 7 


P4i 


4 


3 


hex. perovskite 


CsCuCla 


P6i 


6 



the context of non-symmorphic latices. There, the minimum 
band insulating filling is not two (spinful) electrons per unit 
cell, but a larger integer multiple. Thus, on non-symmorphic 
lattices a fully symmetric insulator with two electrons per unit 
cell immediately implies an exotic ground state, and provides 
guidance to the search for novel phases of matter. We leave 
extensions of these ideas to systems with spin-orbit coupling 
as well as to quasicrystals, to future work. 



TABLE I. Some non-symmorphic groups and their ranks, colloquial 
structure names and representative materials. 



SSL that characterizes a CuB0 3 layer in SCBO (Fig.[3fb)) has 
the space group pAg with a glide symmetry (PigJTJ. In d = 2 
screw axes are impossible so S = 2. Since v = 1 < <S, the 
magnetization plateau must be topologically ordered, or else 
break symmetry as in some proposed candidates.^ 

III. CONCLUSIONS 

We have demonstrated that the textbook reasoning of fill- 
ings at which band insulators may occur must be sharpened in 
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SUPPLEMENTARY MATERIAL 

1. Fourier-Space Approach, 'Many-Body' Crystallography via 
Twist Operators and Electronic Polarization 

In this appendix, we rederive our results (and extend them 
to the two exceptional cases) in a different language, that of 
Fourier-space crystallography. This is a reformulation of crys- 
tallography and space-group classification originally due to 
Bienenstock and Ewa ld and developed more fully by Mer- 
min and collaborators W^MI that is especially well-suited to 
a many-body formulation. As this approach will be unfamil- 
iar to most of our readers, we will summarize the basic ideas 
before discussing the many-body generalization. 

A key idea of the Fourier-space approach to crystallogra- 
phy is to permit all operations that preserve the correlation 
functions of the crystal (i.e., its diffraction pattern) rather than 
its real-space positions — as a consequence, a translation that 
merely shifts the origin of the crystal is equivalent to the iden- 
tity. Building on this, the basic tenets of the Fourier-space 
approach are: 

(i) the action of a group operation on a given Fourier-space 
density (Bragg peak amplitude) is to multiply it by a 
phase factor: 

g : n k -> e 2 ^^n gk (6) 

where $ s (fe) is a linear function of k and the 2ir is in- 
troduced so that $ s (fc) can be chosen rational. 

(ii) a shift of real-space origin acts as a gauge transformation 
on the phase factors: 

X : $ s (fe) -> & g (k) = $ 9 (fc) + X (gk - fc), (7) 

where \ is also a linear function. (Note that this is not to 
be confused with gauge transformations corresponding 
to the conserved U(l) charge.) 

(iii) requiring that the group multiplication rule is faithfully 
represented yields a group compatibility condition on the 
phase factors: 

* flh (fc) = **(M0 + */.(*) (8) 



Note that in order to maintain consistency of notation with 
the original formulation of the Fourier-space approach, we 
define symmetry transformations passively] ' I 13 ! 18 ! in contrast 
to the main text where the transformations are active. The 
classification of different space groups then reduces to the 
classification of all the possible sets of phase factors, up to 
gauge equivalence in the sense above, given a Bravais lattice 
and a compatible point group. Owing to the primacy given 
to the diffraction pattern rather than the real-space structure, 
this formulation of crystallography naturally extends to qua- 
sicrystals. As an aside, we note that in Fourier space crys- 
tallography, it is conventional to speak of 'the' symmetry g 
without worrying about the associated translation, since this 
emerges as the corresponding phase factor $ g (k); really, then, 
the set G ~ {<7|t} in real space is replaced by G ~ {g|3>g} in 
Fourier space, with the caveat that $ 9 is to be understood in 
a gauge-fixed sense — this is equivalent to the statement that 
writing G ~ requires a choice of origin. 

A phase factor $ s (fc) for the action of a group operation g 
on a given Fourier component n k of of the density is gauge 
invariant if and only if the reciprocal lattice vector fc is left 
invariant by g i.e. if gk = k. An essential screw or glide 
has &g(k) ^ (throughout this section, we use '=' to mean 
'modulo integers', since only the fractional part of the phase 
function matters), and manifests itself through the vanishing 
of the Fourier component n k — thus, of the 157 space groups, 
155 are characterized by systematic extinctions in their Bragg 
patterns. The remaining two exceptional cases do not have any 
gauge invariant phase factors; rather, they are characterized by 
gauge invariant phase differences. It is evident that such phase 
differences correspond to an obstruction, not in fixing a gauge 
to render any given screw or glide simple (free of an associated 
fractional translation) but in finding a gauge about which all 
screws/glides are simultaneously simple. Since 'choosing a 
gauge' is Fourier-space language for 'choosing a real-space 
origin', the result follows. Applying symmetry to the Bloch 
theory of energy bands, it was shown by Konig and Mermin 
that the existence of gauge invariant phase differences in space 
groups nos. 24 and 199 leads to a required touching of bands 
in the Brillouin zone, which suggests that this should be true 
in the interacting case as well, and which we confirm below. 

The central result of this appendix — and in some ways the 
most elegant formulation of our results — is that the Fourier 
space approach has a natural interpretation as 'many-body' 
crystallography if we consider the action induced by the sym- 
metry generators on the twist operators. We shall show that in 
this language the phase factors appear in an analogous fash- 
ion to the single particle case, but multiplied by the operator 
Q/N, accounting for the many-body structure; the same no- 
tion of gauge equivalence and the group compatibility con- 
dition will hold. With the usual caveats about the choice of 
N and v, the properties for the ground states demonstrated in 
the main text follow immediately. In addition, the exceptional 
cases appear quite straightforward in this language. This is 
by its nature a rather telegraphic account; we defer a detailed 
discussion to future work. 

We begin by introducing a full set of flux-insertion oper- 
ators. These are defined by first adiabatically turning on a 
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gauge field A = k/N with an integer number of flux quanta 
through each non-contractible loop on the torus; we call this 
operation Fk- Since an integer number of quanta is not ob- 
servable, the final Hamiltonian is equivalent to the initial one, 
after applying a gauge transformation. The necessary trans- 
formation is 



exp 



i 

N 



d d r (fc ■ r)p(r) 



(9) 



The product of the two transformations Tk = UkFk, which 
we will colloquially term a twist operator, returns the Hamil- 
tonian to itself. Furthermore, since the system is assumed to 
be an insulator, the energy cannot change more than an expo- 
nentially small amount. Thus the twist operators {Tk} for re- 
ciprocal lattice vectors fc are operators on the space of ground 
states. Crucially, if the ground state is unique, it must be a si- 
multaneous eigenstate of all of them. This will contradict the 
transformation properties of the {Tk} in a non-symmorphic 
lattice if the filling is not a multiple of the non-symmorphic 
rank S. 

Let us consider a symmetry operator G ~ as before. 

It acts on the gauge transformation operator Uk via conjuga- 
tion: a straightforward calculation yields 

G" 1 i7 fe G = exp|-^ j d d r {gk) ■ rp(r)| e -< fc *~ 1 ' r<s / w 

E^e^W^ (10) 

where we have used the fact that k ■ (</ -1 r) = gk ■ r, and 
identified the phase function 

* g {k) = ^(g- 1 T)-k = ^T-{gk) (11) 

which is manifestly a linear function of k. (Note that when 
gk = k these reduce to the definitions in the main text.) 

The adiabatic turning on of the flux Fk transforms as 
G _1 FfcG = F g k, and so the full flux-threading operator trans- 
forms as 



&- 1 f k 6 = f ak e 2«* B (fc)c/^ 



Furthermore, we have the relations 



k t +k 2 ■ 



(12) 



(13) 



In order to understand the gauge transformation induced by 
a change of origin, consider the transformation r — » r + a 
which accomplishes such a change. We then have 



U k -> U k e 



2Tvi x (k)Q/N 



(14) 



where x(fc) = a ■ k is again a linear function of k. From 
this, it follows that under a gauge transformation we must have 
<& g (fc) -> 3>„(fe), where 



<s>'(k) = <z> g (k) + x (gk)-x(k) 



(15) 



which we obtain by undoing the gauge transformation to re- 
turn to the original gauge, performing the symmetry opera- 
tion, and then redoing the gauge transformation to return to 



the new gauge. Note that ( 14 1 has the same form as (|6j, ex- 
cept that the phase factor is multiplied by the number opera- 
tor. In addition, the phase factor transforms identically under 
gauge transformations — as is evident from comparing ( 15 1 
and (FT). Going further, we can verify that the group compat- 
ibility condition (|8]l is obeyed for a pair of transformations 
{9\r},{h\r'}. 

Thus, the 'many-body' definition of the phase function that 
can be inferred from the action of the symmetry generators on 
the set of twist operators is consistent with that obtained from 
considering the 'classical' action of symmetry on the recip- 
rocal space amplitude of the particle density. Since the phase 
function, the space on which it acts (i.e., the reciprocal lattice), 
the group compatibility condition and the gauge transforma- 
tion all agree, the results of the Fourier-space classification 
of space groups corresponds with a classification of the alge- 
bra of twist operators, with the caveat that the phase factors 
are multiplied by an additional factor reflecting the many- 
body physics. By replacing Q by its eigenvalue N 3 v (valid 
since particle number is a good quantum number), it is evident 
that many-body phase factors differ from the phase factors for 
single -particle Fourier-space crystallography by a factor N 2 v. 
For v = 1 and for an appropriate choice of N 2 that is indivis- 
ible by any of the denominators of the nonzero [i.e., up to 
addition by integers) phase factors, the fact that a phase fac- 
tor or phase difference is nonzero in the single-particle case, 
which flags a non-symmorphic space group, retains the prop- 
erty that it is nonzero in the many-body sense. In this sense, 
the non-symmorphicity of the space group, and attendant con- 
sequences for the phase factors and the symmetry structure, 
carry through to the many-body case. We will demonstrate 
below a few consequences of this correspondence. 

The proof of degeneracy of the 155 non-exceptional non- 
symmorphic groups is quite straightforward in the Fourier- 
space language. Consider the expectation value (4 r |Tfc|4 f ). 
Here, fc is chosen along along the invariant plane/axis of an 
essential glide/screw g, which satisfies g q = 1 (as always, 
q = 1, 2, 3, 4 or 6). Then, $ g (fc) ^ by definition, and since 
we have gk = k the phase is gauge invariant. From this com- 
bined with the fact that g q = 1 and the group compatibility 
condition, we must have q$ g (fc) = 0; thus, the denominator 
of the phase factor must be a divisor of q. A further constraint 
is placed on the phase factor if some power q < q of the phase 
factor is removable, for in this case in some gauge we must 
have q$ g (k) = 0, for otherwise <?' would be essential. If 
we choose q to be the smallest power for which g° is remov- 
able we must have q be a divisor of q (this can be verified by 
explicit computation). Thus, the denominator of 3> ff (fe) must 
be q and the numerator must be relatively prime to this for if 
not, a smaller factor than q would render $ s (fe) = 0. From 
our definition of the non-symmorphic rank of a space group 
operation as the smallest power to which it can be raised to 
make it removable, it is evident that S g = q, and from the 
above discussion we see that the phase factor is of the form 
$ g (fc) = p/Sg with p, S g relatively prime. From this, we 
see that if the ground state preserves the symmetry we have 
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(replacing Q by N 3 v) 

(*|f fc |*> = WG-^fkGW 

= <*|T flfe |*) e «hr*«(*)^" (16) 

Since gk — k and the gauge-invariant phase factor is given by 
$ 9 (fe) = p/<S s , we find that 



(*|T h |*) = (*|T fc |*)e 



i2irN 2 vp/S g 



(17) 



Now if the ground state | ty) is unique, then it is an eigen- 
vector of Tfc, since this operator transforms ground states into 
other ground states. It then follows that the filling must be a 
multiple of S g (by taking TV to be a large number that is rel- 
atively prime to S g ). In other words, if the symmetries are 
preserved, then any sufficiently large system whose size N is 
relatively prime to S g must have degenerate ground states. To 
make the connection to the proof quoted in the main text ap- 
parent, observe that since Tk = UkFk, and Fk, Uk, correspond 
to inserting a flux and performing a gauge transformation, it is 
clear that (*|T fe |*> = (*|?7 fc |*'> = (*|*), and thus the two 
approaches are equivalent (modulo changes due to the active 
versus passive definition used in this appendix.) 

This approach using twist operators can also be used to de- 
fine the non-symmorphic rank arising from the interplay of all 
the symmetries in the group, a particular application of which 
is to rule out a featureless phase for the two exceptional space 
groups not covered by our preceding discussion. To do so, we 
first need to introduce the idea of many-body polarization. A 
sharp characterization of the insulating phase is given by the 
electronic polarization, 4 5 a quantity that is nonzero in the in- 
sulator but ill-defined in the metallic phase. Specifically, in 
this paper we consider a definition of the polarization^ ap- 
propriate to the specific periodic system size N. Recall that 
if the ground state is unique then it must be an eigenstate of 
all the Tfc. Eq. ( 13 i implies that the phase of the eigenvalues 



have a linear dependence on k 



I*) 



(18) 



where we use the superscript to remind us that this definition 



is for a finite system. (We write ( 18 i in units where e = 1; in 



general units, T should be replaced by "P/e throughout.) The 
behavior of the finite-system polarization T>( N > in the ther- 
modynamic limit can be used to define the many-body po- 
larization of the ground state when a Chern number is not 
present. Even in the latter case, although polarization is no 
longer well-defined, "p( N > on a system of a fixed size still 
makes sense; note that as it is a vector, the group operations 
transform "P^ N \ Now, replacing the expectation values in Eq. 



( 16 1 by the exponential of the polarization, and using Eq. ( 1 1 
we find that the combination 



w 



(JV) 



= -uN 2 t ■ gk + V [N) ■ (gk - k) 

= f-l/JVVV + Q-'V^ - V {N) ) ■ k 



(19) 



is an integer. Since this is true for every reciprocal lattice 
vector k, 



g V W + vN 2 T _ V (N) 



where the equivalence is modulo a Bravais lattice vector. Now 
we take ( |20| ) for systems of size N and N + 1, multiply them 
by 2N + 3 and 2N — 1 respectively and subtract. This cancels 
off the factor of N 2 , leaving 



gV + vt = V 



(21) 



where V = (2N + 3)V (N) - (2N - 1)V {N+1) . From this, 
it is clear that 



"P\ V 1 

9\ — )+ T = — + ~ a G 



(22) 



for some Bravais lattice vector ag. Thus, if we define the 
coordinate rp = T 3 jv we find that 



G : r -> gr 



r 



v 



(23) 



(20) 



for some Bravais lattice vector Qg, i.e., r Q is invariant under 
the symmetries up to l/v tYl of a Bravais lattice vector. If v = 
1, then ro can be used as a center for all group operations, and 
hence the lattice must be symmorphic. For a non-symmorphic 
lattice, the minimum value of v for which such a point r 
exists can be defined as the non-symmorphic rank S; this is 
the minimum filling where trivial insulators have a chance to 
exist. 

Intuitively, V is the polarization of a unit cell (relative to the 
origin appearing in the definition of Uk), so Vq = V fv is the 
center of charge of the v electrons in the cell. If there is a sin- 
gle electron per unit cell, this center of charge is well-defined 
(within each unit cell). If the group is non-symmorphic, ap- 
plying the right symmetry will change the charge-center, and 
thus the new wave function is different from the original one. 
When there are more electrons in a unit cell, the center of 
charge becomes less well defined, with v z locations in a grid 
in each unit cell (on account of dividing the polarization, de- 
fined only modulo a lattice vector, by v). Thus, for a large 
enough value of v, it becomes easy to ensure that these charge 
centers are mapped among themselves. 

This argument applies to the two exceptional groups. For 
any one element of these groups, there is a symmetric location 
for the charge center. But the group elements do not have a 
common point of symmetry. We find in fact that 5 = 2 for 
these groups. 

As the remaining 155 non-symmorphic space groups were 
covered by our previous arguments, this concludes the proof 
that all non-symmorphic space groups have S > 1. We ob- 
serve that the phase-space approach motivates two possible 
lines of further inquiry: (i) it extends naturally to higher di- 
mensions, and to quasicrystalsp (ii) it has a natural formula- 
tion in the language of cohomology, 10 and suggests that the 
existence of gapped featureless phases with space-group sym- 
metries may be understood within this language. Similar tech- 
nology has been brought to bear, with fruitful consequences, 
for situations with local symmetriesP^l 
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2. Examples in the Main Text 

a. Shastry •Sutherland Lattice 

The SSL is generated by the primitive vectors ai = 
(b + c)x,a 2 = (b + c)y and the four-site basis 
{l(x-y),l(-x + y),§(x + y),-l(x-y)} with b + 
c, and corresponds to 2D space group pAg. The non- 
symmorphic operation is a glide comprised of a reflection 
about a plane parallel to a\ and passing through the origin 
in the specified coordinate system, followed by a translation 
by ai/2. Since the two primitive vectors are orthogonal, it 
follows that this cannot be the projection of a lattice vector on 
the glide plane and therefore the glide is essential, leading to 
S = 2 (in two dimensions, this is the only nontrivial value 
of S). Since there are four sites per unit cell, a spin-| model 
will have saturation magnetization of 2 per unit cell; half this 
value corresponds to a filling of v = 1, and thus the half- 
magnetization plateau must either break symmetry or exhibit 
topological order. 



b. Hexagonal Close-Packing 

The hep lattice is generated by the primitive vectors a\ = 
ax,a 2 = a (^—^x + ^y^J ,a 3 = cz and two-site ba- 
sis {0, |ai + \a 2 + and corresponds to space group 
P6s/mmc (no. 194). The crystal possesses a screw axis 
involving a sixfold rotation about a 3 centered at the point 
(2ai + a 2 )/3, followed by a translation by \a 3 . Applying 
this twice yields the product of a pure translation (by a 3 ) and 
a three-fold rotation, and so the non-symmorphic rank of this 
operation is 2. The only other non-symmorphic element is a 
glide reflection, also of rank 2, from which we conclude that 
the rank of the group is S = 2. Our result then shows that 
in the absence of broken lattice symmetries, a noninteracting 
band insulator (of spinful electrons) must be metallic unless 
the filling is a multiple of 4. Furthermore, if interactions are 
included and open a gap, except at such fillings the result- 



ing insulating ground state cannot be featureless, and so must 
break symmetry or exhibit topological order. 

c. Diamond 

The diamond lattice is characterized by the non- 
symmorphic space group Fdim. In coordinates in which 
the fee Bravais lattice primitive vectors are ai = |(y + 
z), a 2 = §(£ + x), a 3 = |(y + z) and the two-site basis 
is {0, j(x + y + z) }, the crystal has a 4i screw axis consist- 
ing of a fourfold rotation about the z axis centered at x/A, 
followed by a translation through £/4. Applying this twice 
yields a two-fold screw that is removable, as can be readily di- 
agnosed by the fact that lattice vectors ai,a 2 have projections 
of z/2 onto the axis. As all the remaining non-symmorphic 
operations are glides which also have rank 2, S = 2. Now, 
take a model of local spin-^ moments on the diamond lat- 
tice, described by a Hamiltonian with U(l) spin symmetry and 
hence conserved magnetization along some axis, which we la- 
bel S z . Following our prescription, we identify the charge at 
position r as Q r — S z + \. From this, we find that the total 

charge is Q = £ r Q r - S* + iV 3 , where S z ot = £ r S z and 
we have used the fact that there are two spins in each unit cell. 
This corresponds to a filling of v = s z + 1 where s z is the 
magnetization per unit cell. A gapped, featureless, and topo- 
logically trivial paramagnetic ground state (the spin analog of 
a featureless insulator) is ruled out unless v is a multiple of 
2. In particular, if in addition to the U(l) symmetry the Z2 
symmetry S z — > — S z is preserved (which includes the fully 
SU(2)-symmetric limit) the total magnetization must vanish 
and so v = 1. In this case, any gapped ground state must 
either break lattice symmetry, or form a topologically ordered 
spin liquid. This suggests a route to a three-dimensional quan- 
tum spin liquid on the diamond lattice by frustrating the ten- 
dency of spins to magnetically order. We also observe that the 
pyrochlore lattice has the same space group as diamond, but 
four sites in each unit cell; for a spin-| model on this lattice 
with vanishing total magnetization, we find v — 2 and there- 
fore a gapped, featureless quantum paramagnet is an allowed 
ground state. 
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